Anomalous In-Plane Anisotropy of the Onset of Superconductivity in (TMTSF) 2 C10 4 
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We report the magnetic field-amplitude and field-angle dependence of the superconducting onset temperature 
jonset Q f me or g an i c superconductor (TMTSF) 2 C104 in magnetic fields H accurately aligned to the conductive 
ab' plane. We revealed that the rapid increase of the onset fields at low temperatures occurs both for H || V and 
H || a, irrespective of the carrier confinement. Moreover, in the vicinity of the Pauli limiting field, we report a 
shift of a principal axis of the in-plane field-angle dependence of T° nxt away from the b' axis. This feature may 
be related to an occurrence of Fulde-Ferrell-Larkin-Ovchinnikov phases. 
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Since the discovery of the organic superconductors 
(TMTSF) 2 X (where TMTSF stands for tetramethyl- 
tetraselena-fulvalene, X = C10 4 , PF 6 , etc.) 0], their 
superconductivity has been studied with much attention. 
Because of the strong anisotropy in the electrical conductivity 
of these materials yfl, they provide excellent opportunities 
to study the properties of quasi-one-dimensional (Q1D) su- 
perconductors. One of the most important and controversial 
issues on the superconductivity of these materials is their 
superconducting (SC) pairing symmetry fl]. In this Letter, 
we provide experimental results that contain new crucial 
clues in understanding the SC symmetry of (TMTSF) 2 C104. 

It has been suggested that the superconductivity of 
(TMTSF) 2 X is unconventional with line(s) of node on its SC 
gap, indicated through the NMR relaxation time 0] and the 
impurity concentration dependence of the transition tempera- 
ture T c JfJ. However, its SC symmetry is still controversial 
as we will review below. One key feature of the SC symme- 
try is their unusually-high upper critical fields H c o{T). Lee 
et al. reported that H c2 (T) of (TMTSF) 2 PF 6 determined 
from resistivity diverges as temperature decreases and H C 2(T) 
reaches up to 80 kOe at the lowest temperatures when mag- 
netic fields H are applied parallel to the b' axis (perpendic- 
ular to the most conductive a axis in the ab plane). In this 
field direction, carriers are confined in the ab plane due to 
the field-induced dimensional crossover (FIDC) as the field 
increases flE!]. This confinement can be interpreted as a con- 
sequence of cyclotron motion of carriers on the Fermi surface 
(FS) 0]. The FIDC suppresses the orbital pair-breaking ef- 
fect and may allow the superconductivity to survive in higher 
fields. Interestingly, 80 kOe for H c2 II V far exceeds the 
so-called Pauli-Clogston limit Hp [10], which fulfills a rela- 
tion Hp/T c = 18.4 kOe/K for an isotropic gap, where sin- 
glet Cooper pairs are unstable because unpaired carriers have 
a lower energy due to the Zeeman effect. In the case of 
Ref. |01, Hp was estimated to be 20 kOe. Similar results 
have been obtained in (TMTSF)2C104 by resistivity and mag- 
netic torque measurements ll ill . One interpretation attributes 
this survival of superconductivity above Hp to a spin-triplet 
state 11121 1 1311 . On the other hand, it has been pointed out 



that, in Q1D superconductors, even singlet superconductivity 
can be stable far above Hp by forming a spatially-modulating 
SC state |[l4l [l5l 161. which is equivalent to the so-called 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state B17L llal . In 
2002, Lee et al. 

reported the absence of a change in 
the 77 Se Knight shift of (TMTSF) 2 PF 6 at T c under pressure, 
in favor of a triplet scenario. However, recently Shinagawa 
et al. ll20ll observed a clear change of the 77 Se Knight shift 
of (TMTSF)2C104 at T c in lower fields. This finding mo- 
tivated us to reexamine the possibility of singlet pairing in 
(TMTSF) 2 X. 

To resolve this puzzle, we are interested in the superconduc- 
tivity in H || a and its in-plane anisotropy. Although not much 
attention has been paid to the superconductivity for H \\ a so 
far, data for H c2 (T) || a of (TMTSF) 2 PF 6 looks quite in- 
teresting: It has a steep slope near H — 0, but saturates when 
it reaches Hp probably due to the Pauli effect, and it slightly 
increases again below 0.3 K. However, for (TMTSF)->C104 
H C 2(T) || a was reported only above 0.5 K 112 ill . The in-plane 
anisotropy of H c2 of (TMTSF) 2 C104 was also reported but 



only at 1.03 K EM, where H c2 (T) is far below Hp. 

In the present study, we revealed the rapid increase of the 
onset fields not only for H \\ b', where the electronic state 
becomes essentially 2D due to the FIDC, but also for H || a, 
where the electronic state remains anisotropic 3D. We also 
observed new features of the in-plane anisotropy developing 
above 20 kOe, which provide a crucial step to understand the 
origins of the enhancement of H c2 , in terms of FFLO states. 

We used single crystals of (TMTSF) 2 C104 grown by an 
electro-crystallization technique, with dimensions of approx- 
imately 2.0 x 0.2 x 0.1 mm 3 . We report here the results of 
the sample with the highest T c among up to 10 samples. We 
note that we obtained similar results in another sample. The 
resistance along the c* axis R c . was measured using an ac four- 
probe method. (The direction of the c* axis is perpendicular to 
the ab' plane and is the least conductive direction.) The mea- 
surements were performed with a dilution refrigerator down to 
80 mK. Temperature was measured using a Ru0 2 resistance 
thermometer with magnetoresistance correction. The anion 
ordering temperature of (TMTSF) 2 C104 is 24 K. Therefore, 
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FIG. 1: (color online) (a) Temperature dependence of R c » in zero 
field. The directions of the orthogonal crystalline axes are illus- 
trated in the inset, (b) Conductance difference Act = R^iH^ = 
0) - R~}(H C - > 0) under an in-plane 50-kOe field applied parallel 
to the b' axis. Resistances for H c * = kOe and H c , = 0.5 kOe 
are plotted against the right vertical axis, (c) Temperature depen- 
dence of Atr with the in-plane field applied parallel to the a axis and 
H c - = 1.0 kOe. Some data are shifted vertically for clarity. 



in the temperature interval between 25 K and 22 K, a cooling 
rate as slow as 2 mK/min was chosen to ensure that all anions 
are ordered and the whole sample is in the "relaxed state". 

Magnetic fields are applied using the "Vector Magnet" sys- 
tem M22I1 . with which we can control the field direction without 
mechanical heatings. The directions of the orthogonal crys- 
talline axes (the a, b', and c* axes) of the sample were deter- 
mined from the anisotropy of H C 2 at 0. 1 K. The accuracy of 
field alignment with respect to the ab' plane and of the a axis 
within the ab' plane are both better than 0. 1 degree. We also 
determined the directions of the triclinic crystalline axes (the b 
and c axes) from angular magnetoresistance oscillations. The 
details of these procedures will be presented elsewhere. Here- 
after, we denote the azimuthal angle within the ab' plane as <p 
which is measured from the a axis. We defined <p so that the b 
axis lies in the quadrant 0° < <p < 90° as indicated in Fig. [3] 

We first present R C -(T) in zero field in Fig. \V[a). Although 
R c - of this sample started to drop at as high as 1.45 K and 
reached zero at 1.30 K, R c - increases again below 0.8 K. This 
increase, which is almost independent of magnetic fields, is 
probably attributed to small cracks in the sample. The data 
of RAT) for H || b' at 50 kOe is presented in Fig. Htb). 
We observed a decrease of R c > below 0.2 K, consistent with 
a previous report 12311 . In order to confirm that such a de- 
crease is due to a superconducting contribution, we mea- 
sured R C -(T) after adding a small out-of -plane component 
H c > = 0.5-1.0 kOe to the magnetic field. If this decrease is 
due to the superconductivity, H c * should suppress the super- 
conductivity and eliminate the decrease of R C -(T). As plotted 
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FIG. 2: (color online) Magnetic fields vs. r° nset of (TMTSF) 2 C10 4 
for H || a (filled squares) and H || b' (filled circles). The arrow indi- 
cates the estimated value of Hp = 26.7 kOe. The phase diagram for 
H || c* is shown in the inset. The broken lines indicate the calculated 
initial slopes using transfer integrals explained in the text. 



in Fig. [TJb), the decrease was indeed eliminated by adding 
H c >. Therefore, it is confirmed that the decrease of R C >{T) 
is a contribution of the superconductivity. We used the fol- 
lowing procedures to define the onset temperature of super- 
conductivity r° nset : We evaluated the conductance difference 
Act = R- l (H c . = 0) - R-}(H C , > 0) and defined r° nset as 
the temperature at which Acr(T) exhibits a sharp increase, as 
marked by the small arrow in Fig.[TJb). This definition char- 
acterizes the very onset of superconductivity. We note that 
this anomaly in Acr(T) is not due to the normal state magne- 
toresistance, because it is unlikely that an abrupt change in 
the difference between R C (H C * = 0) and R C -(H C * > 0) oc- 
curs at a certain temperature. The definition has the advan- 
tage that T° nset is not affected by the extrinsic small increase 
of R c * because it is cancelled in the subtraction. For H || c*, 
T° nset (H) was determined similarly from the conductance dif- 
ference Acr(H) = R-}(H C , = H) - R-}(H C , = H + AH). 

The phase diagrams for H || a, H \\ b' , and H \\ c* are pre- 
sented in Fig. [2] In the vicinity of H — 0, linear temperature 
dependences of the curves were observed for all field direc- 
tions, which can be analyzed in a GL model for a clean type- 
II superconductor. Within a GL theory with a tight-binding 
model, the slope dH C 2{T)/dT at T C (H = 0) is related to the 
transfer integral t of each direction J24{]. By taking into ac- 
count the k z dependence and the nodes of the gap over the FS, 
we obtain t a = 1200 K, t b , = 310 K, and t c . = 7.0 K from the 
initial slopes indicated by the broken lines in Fig. [2] These val- 
ues agree favorably with realistic band parameters J^t]. From 
these analysis, it is clear that H C 2(T) is governed by the orbital 
limitation at low fields in all three directions. 

In higher fields, the behavior of these curves is qualitatively 
different. The curve for H \\ V keeps a linear temperature de- 
pendence up to 35 kOe and starts to exhibit a rapid upturn 
in higher fields. This behavior is consistent with the "initial 
cool" curve in Ref. Ill 111 . For H || a, the curve apparently 
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FIG. 3: (color online) Polar plots of T° nset (4>). The points for |0| > 90° 
are from the same data as those for \<p\ < 90°. They are plotted in 
order to symmetrize the figure. Error bars are omitted for the sake 
of clarity. Purple crosses indicate that no SC anomaly in Acr(T) was 
observed above 80 mK. The directions of some crystalline principal 
axes are illustrated in the upper-left graph. Solid red lines indicate 
the direction of the new principal axis X. 



shows limiting behavior; this is consistent with the Pauli- 
limiting behavior with the estimated value of Hp = 26.7 kOe. 
Interestingly, in higher fields, a small kink of Act remains vis- 
ible up to 50 kOe also for H \\ a, as shown in Fig. [He). Con- 
sequently, the onset curve diverges for H || a at low temper- 
atures, too. This is, to our knowledge, the first report of the 
low-temperature high-field phase diagram of (TMTSF^CICU 
for H || a. We note that we obtained similar phase diagrams 
for another sample. It is interesting that all three onset curves 
in Fig. |2]look similar to those for (TMTSF) 2 PF 6 Q. 

Next, we focus on how T° nset changes when magnetic fields 
are rotated in the ab' plane. The data are displayed in Fig. [3] 
using polar plots of T° nset (<p), where the direction of each point 
seen from the origin corresponds to the field direction and the 
distance from the origin corresponds to T° nset . At low fields, 
T° nsel ((p) exhibits a sharp cusp at (p = 0° (H || a) and a broad 
minimum around </> = +90° {H \\ b'). These low-field results 
are consistent with H C 2(<P) reported by Murata et al. 12111 . al- 
though the sharp peak at (f> — 0° cannot be explained in an 
anisotropic 3D GL theory 12511 , The chain-like crystal struc- 
ture may play an important role in generating the sharp peak. 

As the field increases above 20 kOe, another anomaly i, e. 
dips of r c onset (0), emerges at |^| = dip = 17 ± 1°. We note that 
R C '(T) in the normal state exhibits non-metallic temperature 
dependence for \(p\ > (pj,o-2n- 19 + 1° above 20 kOe, signal- 
ing the onset of the FIDC SIH- Because ^ip — 03D-2D, we 
infer that these dips are related to the FIDC. When the dimen- 
sionality of the electronic system is lowered, superconductiv- 
ity in in-plane magnetic fields becomes more stable because 
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FIG. 4: Field dependence of the relative difference between 
r° nset (+45°) and r™ set (-45°). Appearance of the new principal axis 
X results in finite values of this quantity. The arrow indicates the 
estimated value of Hp = 26.7 kOe. 



the orbital pair-breaking effect is suppressed. Thus T° met ((p) 
should be enhanced for \(f>\ > 03d-2d> resulting in a minimum 
of r c onset (0) around 3 d-2d- 

The third and probably the most important anomaly is, that 
in magnetic fields above 30 kOe, the V axis is no longer a 
symmetry axis of r° nset (</>) and a new principal axis X appears 

around (p 70° as indicated by the solid red lines in Fig. [3] 

Moreover, behavior of r° nset (</>) around X, a principal axis at 
high fields, and b', a principal axis at low fields, is qualita- 
tively different: At high fields T° nset {(p) is enhanced around 
X, while at low fields T° nset (<p) exhibits a broad local mini- 
mum around the b' axis. In addition, this X axis tends to rotate 

from (f> 70° toward the b' axis as the field increases. At the 

largest field available in this study, the deviation of X from the 
V axis is reduced to about 10°. We checked that this change 
of symmetry is not due to a misalignment of the magnetic 
fields. In Fig. |4l we plotted the relative difference between 
r c onset (+45°) and r c onset (-45°) against the field strength. This 
quantity represents the asymmetry with respect to the b' axis, 
thus the appearance of X results in finite values. It is evident 
that the asymmetry, ;. e. X, is absent in lower fields and then 
starts to develop around Hp. Therefore, the appearance of X 
cannot be attributed to conventional origins like an anisotropy 
of the Fermi velocity, because variation of r° nset (</>) from such 
origins should develop from H — 0. 

We now discuss the origin of the new principal axis X, 
which also indicates the direction of the field in which 
T° nset {(p) is enhanced. Its appearance should be related to 
the Pauli pair-breaking effect, because X appears at nearly 
Hp. In the case of singlet pairing, the appearance of X is at- 
tributable to the formation of an FFLO state Il7lll8ll . in which 
the Cooper pairs have a finite wavevector ^fflo- In a Q1D su- 
perconductor, the stability of this state is greatly enhanced by 
the nesting properties of its FS J26ll and that ^fflo essentially 
matches the nesting vector between the spin-up and the spin- 
down FSs, which should be nearly parallel to the a axis and 
should be independent of the field direction. 

For H || b', it has been discussed using orthorhombic band 
structures that an FFLO state with </fflo II a becomes stable 



with a help of the FIDC H [3 [I3l. Although we are not 
aware of a publication on the in-plane field-angle variation of 
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this FFLO state taking into account the realistic triclinic band 
structure of (TMTSF^CIO^ we expect that for \<p\ > </>3d-2d, 
where the FIDC takes place, this FFLO state is still stable. 
However, the direction of ^fflo, matching with the nesting 
vector, should be slightly tilted from the a axis because of the 
triclinic FS of (TMTSF^CICU. In addition, ^fflo may vary 
with increasing the field because the separation between the 
spin-up and the spin-down FSs depends on \H\. Within this 
scenario, one possible explanation of X, which is the field- 
dependent special direction in the range \(p\ > 03D-2D, is that 
X is perpendicular to ^fflo and thus H \\ X corresponds to 
H ± <7fflo ■ Because the direction of ^fflo is expected to 
depend on the field strength as we explained, X (_L ^fflo) 
may also rotate in increasing field, which is consistent with 
our experimental results. 

While for \<p\ < 03D-2D, namely near H || a, the nature of 
superconductivity may differ from that near H || b' because of 
the absence of the FIDC. Despite the absence of the FIDC, the 
orbital-limiting field is much larger than Hp near H || a at low 
temperatures, which is evident from the steep slope of H C 2{T) 
at H = in Fig. [2] An FFLO state in a Q1D system in fields 
parallel to the most conductive axis has been proposed in a 
study of doped two-leg ladder cuprates using a t-J model 12711 . 
We infer that a similar FFLO state might be stable near H \\ a, 
although a theory adapted to (TMTSF^QO/t, a coupled chain 
system, needs to be developed. The recent NMR study, which 
reported that the density of states at the Fermi level recovers to 
the normal state value in the SC phase above 20 kOe for both 
H || a and H \\ b' ll20Tl . would support these FFLO scenarios. 

On the other hand, if (TMTSF)2C104 is a triplet supercon- 
ductor, polarized Cooper pair spins may cause an anisotropy 
of T° nset (<fi). Assuming that the spins of the Cooper pairs are 
fixed to one direction, superconductivity is not affected by a 
Pauli effect when the field is exactly parallel to the spins, while 
it is suppressed for the other field directions. In this case, how- 
ever, it seems difficult to explain the rotation of X. 

In summary, we have studied the in-plane anisotropy of the 
superconducting onset temperature r° nset of (TMTSF)2C104. 
We observed that !T° nset remains finite up to 50 kOe in mag- 
netic fields parallel to the a axis, as well as for H || V . We 
suggest that the field-induced dimensional crossover plays an 
important role for the enhancement of T° aset (<f>) when the field 
is tilted more than 17° from the a axis. In addition, we noticed 
that one of the principal axes for superconductivity, which 
points along b' at low fields, shifts away from this direction 
around 30 kOe but evolves back toward the b' axis at higher 
fields. The survival of superconductivity far above Hp and the 
unusual in-plane anisotropy observed in the high field regime 
suggest the stabilization of modulated superconducting phases 
when high fields are aligned to the ab' plane, in favor of a spin- 
singlet scenario. We speculate that two kinds of FFLO states 
are realized in this compound: the one predicted by Dupuis 
et al. 115I1 near H \\ V and the one related to the prediction by 
Roux et al for H \\ a, separated by the dips of T° nsel ((f>) 
around ~ +17°. We believe that theoretical studies taking 
into account the triclinic band structure are desirable to under- 



stand our results and reveal the SC symmetry of (TMTSF^X. 
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